Schwinger Pair Production at Finite Temperature in QED 
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We use the evolution operator method to find the Schwinger pair-production rate at finite tem- 
perature in scalar and spinor QED by counting the vacuum production, the induced production 
and the stimulated annihilation from the initial ensemble. It is shown that the pair-production 
- - - rate for each state is factorized into the mean number at zero temperature and the initial thermal 

' distribution for bosons and fermions. 
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Vacuum polarization and pair production have been issues of continuous concern since the early works by Sauter, 
Heisenberg and Euler, and Weisskopf and then by Schwinger 0] (for a review and references, see Ref. Q). The task 
'"pi ■ of directly computing, without relying on the electromagnetic duahty, the effective action in electric field backgrounds, 
' [ however, has been a challenging problem due to the vacuum instability. Dunne and Hall used the resolvent method to 
directly find the effective action in time-dependent electric fields [1] . In the previous paper [5| , employing the evolution 
' operator method, we found the exact one-loop effective actions of scalar and spinor QED at zero temperature in a 
I— 'i constant or a pulsed electric field of Sauter-type, which satisfy the exact relation 2Im£off = ± J2n ln(l ± A/",!) (with + 
for scalar and — for spinor) between the imaginary part of the effective Lagrangian density Ccs and the mean number 
of created pairs 77n at state n. Even finding the pair production rate by time-dependent or spatially localized electric 
fields is methodologically nontrivial, which has recently been intensively studied [1, 0, Hi • 
I To calculate the effective action and thereby Schwinger pair production at finite temperature is another challenging 
' problem in QED. In a constant pure magnetic field the QED effective action was studied at finite temperature [9| 
' and at finite temperature and density However, the presence of an additional electric field raised Schwinger 

, pair production at debate depending on the formalism employed. Some calculations in thermal field theory reported 
T-H ' that the effective action in both a constant magnetic and electric field had an imaginary part having dependence on 
00 , temperature [HI, IH, ■ However, in the real-time formalism, no imaginary part was found in the QED effective action 
' in the presence of both magnetic and electric field [l3| • The recent calculation of effective action in the imaginary-time 
^ ' formalism shows an imaginary part only at two-loop but not at one-loop in both a constant electric and magnetic field 
• p^ . In nonequilibrium quantum field theory of scalar QED, a calculation in the real-time formalism shows thermal 
r^J I enhancement of pair production [16j]. 

^ . In this paper, using the evolution operator method, we find the pair-production rate at finite temperature in time- 
■ " " ' dependent electric fields both in scalar and spinor QED. The evolution operator, unitarily transforming the particle 
and antiparticle operators from the ingoing vacuum to the outgoing vacuum, carries all the information of quantum 
evolution. In fact, the evolution operator is completely determined by the Bogoliubov coefficients. The advantage of 
the evolution operator is the readiness to calculate the probability for transitions among multiparticle states. This 
allows us to compute the mean number of created pairs at finite temperature in scalar and spinor QED by counting 
the pairs from the vacuum and the induced production and the stimulated annihilation in a thermal ensemble of 
bosons and fermions. We find that the mean number of created pairs is factorized into the mean number of created 
pairs at zero temperature and the initial thermal distribution for bosons and fermions. 

The organization of this paper is as follows. In Sec. II, rewriting the Bogoliubov transformation as a unitary 
transformation by the evolution operator, we find the mean number of created pairs at zero temperature in scalar 
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and spinor QED. In Sec. Ill, we calculate the mean number of created pairs at finite temperature and apply it to the 
Sauter-type electric field. 

II. EVOLUTION OPERATOR AND PAIR PRODUCTION AT T = 

We first consider scalar QED for spinless charged bosons under an external electric field with the gauge field A^. 
The electric field is assumed to be acting on for a finite period of time so that the ingoing and the outgoing vacua 
are well-defined. Thus, at tin = —00, before the external electric field being turned on, the scalar field is free and the 
Hamiltonian takes the usual form 



= j (^;;^^k,i„Afk,in = j ^^;^Wk,inK,inak,i„ + 0^,i„&ka, 



(1) 



where Wk.in is the initial frequency at momentum k. Here, the gauge is chosen = 0, so that the ingoing vacuum 
|0;tin) is the Minkowski vacuum |0)m, annihilated by ak(tin) and fok(^in) for each momentum k. Similarly, the 
outgoing vacuum at tout = 00 is defined by ak(^out) and t'k(iout)- These operators are related through the Bogoliubov 
transformations 

Qk.out = Mkak.in + i'k^kjn' ^k.out = Mk&k.in + ^^kaik.in' (2) 

where |^kp - It'kP = 1- 

To express the outgoing vacuum as multiparticle states of the ingoing vacuum, we rewrite the Bogoliubov transfor- 
mations (m as a unitary transformation 5| 

ak,out(A) = ;7k(A)ak,i„(0)[/^(A), 6k.out(A) = (7k(A)6k,i„(0)(7i[(A). (3) 
Here, the evolution operator C/k is factorized into the overall phase factor and the two-mode squeeze operator as 

where 

e^'^-^^, a-^, 7k = -21n(^k). (5) 

Mk Mk 

From the charge neutrality of the vacuum, equal numbers of particles and antiparticles are produced at zero 
temperature and they carry the opposite momenta due to the momentum conservation. The multiparticle state of 
n-pairs can be concisely denoted as [nk, <} = l?^k; "-k; The probability for n-pairs with momentum k to be created 
from the vacuum is 

P„(k) = |(nk,out|0,in}p = |(nk,in|;7i[|0,in)|2 = e-'-'iekl'". (6) 

Note that Pa = e<^ and Pi = eT"-|a|2 so that F„ = Po(A/^o)" and ErT=o ^" = 1 ^o^' ^ach k. Thus, at zero 
temperature, the mean number of pairs created from the vacuum for each momentum per unit volume is 

00 

Ark(r = o) = ^nP„(k)^|^kr (7) 

Next, in spinor QED, before the interaction of an external electric field, the spinor field is free without the gauge 
potential {A^^ = 0), and has the Hamiltonian given by 

-f^i^nMN^M = H y (^^«an(^Un^«an + 4,1.1^^",™), (8) 

where 6„^in and d„.in are particle and antiparticle operators in the ingoing vacuum. After the interaction of the 
electric field, the ingoing vacuum evolves to the outgoing vacuum, whose particle and antiparticle operator are &„,out 
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and dn,out- The Bogoliubov transformations between the ingoing and the outgoing particle and antiparticle operators, 
hn,dn, are similarly given by 

6n,out = Ain^n.in + i'^n'A.m^ rfn,out = fJ^ndnM - (9) 

where |/inP + l^^nP = 1 and n = (k, cr) with a = ±1/2. The Bogoliubov transformation can be also written as a 
unitary transformation (loj 

^n.out — f^n^n, in ^7^5 ^n,out — UndnSriU^^ (-10) 

where 

Here, the three parameters £,n,ln and On are determined by the Bogoliubov coefhcients as 

?n = -*^, 7n = -21n(|^„|), 62'"^"=^. (12) 

Note that the pair production on spinor QED is restricted to only one pair of particle and antiparticle for a given 
quantum number n due to the Pauli exclusion principle. Thus, the mean number of pairs created from the vacuum 
for each state n at zero temperature is calculated as 

jr^{T^Q) = |(l„,out|0,in)|2 = |^„|2. (13) 



III. PAIR PRODUCTION AT T / 

We now calculate the mean number of pairs at finite temperature from the probability for each transition, as in 
the case of zero temperature. As there is no mode-mixing, we separately calculate the mean number of created 
pairs for each mode. For an initial thermal ensemble at /? = 1/kT, which might not be charge neutral, the mean 
number of produced pairs consists of the vacuum pair production, the induced pair production, and the stimulated 
pair annihilation as shown in Fig. 1: 



.n„>0 p„>n„>0,g„>m„>0 



(14) 



where Z„ is the partition function for the initial ensemble and Pn„,m„^p„.q„ — |(Pn, in|?7n|n„, m„, in)p is the 
transition probability from |n„,m„,in) to |p„,g„,in). Using |(p„, 9„, in|C/„|n„, m„, in) p = |(p„, g„, in|J7t |n„, m„, in) p, 
which implies Pp^,g^_>„^_,„^ = Pn^_rnn^pn.q„, the mean number of pairs created at finite temperature, Eq. (|14p . can 
be written as 



Ar„(T) 



2Zn 



E 



(n„,m„, in|e ''-'^"■'"[/„|p„, g„, in) (p„, g„, in|(iV„anC^rl - U^Nn 

,in ) l^n 1 ■, HI/ 



2Z, 



- (n„, m„, in|[/^e"''-"^"-'°|p„, g„, in) (p„, g„, in|(7V„anC^„ - UnNn:in)\nn, ?ti„, in) 
-TriUnAnNnMUl i„ - iV„,i„)e-'3^".'», 



(15) 



where Nn,\n — ajL^n'^kan + ^k.in^k.in for scalar particles and Nn,\n — ^I.in^n.in + AiM-i^nM for spinor particles and 

Hn-m = i^n^inNnAn for both. 

From the Bogoliubov transformations ([3]) and (|10[) . we have Nn,ont — UnNnAnU}^. Thus, the mean number of pairs 
created at finite temperature can be written concisely as 



A^n(T) - — Tr(iV„,out-^n,i„)e-'3^"-'". 



(16) 
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m„>0 
(a) 



m„ + n„^Q, P „^ q „- m „ 
(b) 



p„>n„>0, Q„>m„>0 
(c) 



FIG. 1: (a) the vacuum production, (b) the induced production, and (c) the stimulated annihilation 
First, for scalar particles, from Eq. ([2]), 

^k.out - ^k,in = 2|iykp(aj^ in«k,in + &k in^k,in + 1) + 2/Xki^kak,in&k,in + 2A<k'^kak.in^k in' 



(17) 



and Zk = e'^'^''■■"(2sinh/3wk,in/2)"^ the mean number of scalar pairs created in momentum k at finite temperature 
is given by 



Ar^^(r)^^ (n,m|(n + m + l)e-("+™)'5"'<^'"|n,TO) = |i.k|'coth^^^ 



A7'^P(T) = --^ (n,m|(n + m- l)e-("+")''"-'"|n,m) = |i/„|2tanh 



(18) 



Second, for spinor particles, from Eq. ([S]), 

A^«,out - ^n,in = -2|i^„|^(6,\ ijj&„,in + d\ -^^dn,in - 1) + 2l^„I/„6„,ind„,in + 2i^* 1^*6]^ j^dj, j^^, (19) 

and Z„ = (l + e-'^'^"'-)2, we find the mean number of spinor pairs created in state n at finite temperature 
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(20) 
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FIG. 2: The mean number of scalar(dashed line) and spinor(solid line) pairs created, A^, in k = state for a pulsed Sauter- 
type electric field of duration t= 10: (a) dependence on temperature T with E= 0.5, and (b) dependence on electric field 
strength E with T= 0.5 for the upper and the lower lines, and with T= for the middle lines, where dimensionless parameters 
ksT /m(? ~+T, qEh/m^c^ — >E, and mTC? /h — »t are used. 



As an interesting model for discussions, we consider the Sauter-type electric field £{t) = £^sech^(t/T) with the 
gauge choice, Az{t) — —Et{1 + tanh(i/r)), which allows the exact solution leading to the exact mean number of pairs 



created, |fkp7 at zero temperature both in scalar and spinor QED 
of the WKB or worldline instanton method has been developed 



IGI . [20I . I21I . Also the approximation scheme 



8|. Because the Sauter-type electric field acts 
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effectively for a finite period of time r, the ingoing thermal states are stable and well defined as required in our 
formalism. Fig. 2 shows the mean number of created pairs Eqs. (|18l20p at finite temperature with the electric field 
E, the duration r and temperature ksT — 1/ (3 scaled in terms of the critical strength, the Compton time and the 
electron mass, respectively. The production of scalar pairs is thermally enhanced, while the production of fermion 
pairs is thermally suppressed as expected by the Pauli blocking, which is consistent with the calculation by density 
matrix method [22] ■ 

The pure thermal effect on the mean number of created pairs, which is /!!j^\s_{T) — 77]i{T) — A7'k(0), is given by 
AJ>7k{T) = ±2|t'kp/if (T) with f^(T) being the Bose-Einstein or Fermi-Dirac distribution. For ksT < cjk.in, the 
distribution f^{T) approximately equals to the Boltzmann factor /k « g-N/m^+k^/fc^T^ xhus, the mean number of 
created pairs at finite temperature, Eqs. p8l20p . is reduced to the zero temperature result, Eas. (|7ll3p . for a temperature 
much lower than the rest mass. 

Our result Eq. (j20p could be compared with the calculation in imaginary-time formalism , where the thermal effect 
appears only at two-loop because thermal one- loop fluctuations are on-shell. On the other hand, Eq. ((20| . an off-shell 
calculation, is equivalent to the thermal loop times the vacuum one-loop, in fact, part of two-loops. One interesting 

comment to be pointed out is that the momentum integral over the distribution function in AA/'k(T) = ±2|fkp/k C^) 
leads to the factor in Ref. , though the momentum integral of the distribution is intertwined with the vacuum 
pair-production rate. To show rigorously the connection between our result and Ref. requires calculating the 
effective action at finite temperature along the line of Ref.j5j, which will be addressed in the future. 



IV. CONCLUSION 



In this paper, using the evolution operator method, we found that the mean number of created pairs in state n at 
finite temperature is given by 

7^n{T) = -^Tr(iV„,out-A^n,i„)e-'3^"-. (21) 

For scalar and spinor QED in external electric fields, the total mean number density of created pairs at finite tem- 
perature is given by 

jr\T) = y'i!^Ark(0)coth^ (for scalar), (22) 
^"(T) = Y.j ^A^f(0)tanh^ (for spinor), (23) 

where A7'^'^(0) and '7T'^{G) are the mean number of created pairs at zero temperature for scalar and spinor QED, 
respectively. 



Acknowledgments 

S. p. K. would like to appreciate the hospitality of Hanyang University, and H. K. L. and Y. Y. would like to 
appreciate the hospitality of Kunsan National University. The work of S. P. K. was supported by the Korea Research 
Foundation Grant funded by the Korean Government (MOEHRD) (KRF-2007-C00167) and the work of H. K. L. 
was supported by the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korea government 
(MOST) (No. ROl-2006-000-10651-0). 



[1] F. Sauter, Z. Phys. 69, 742 (1931); W. Heisenberg and H. Euler, Z. Physik 98, 714 (1936); V. Weisskopf, K. Dan. Vidensk. 

Selsk. Mat. Fys. Medd. XIV, No. 6 (1936). 
[2] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[3] G. V. Dunne, "Heisenberg-Euler Effective Lagrangians: Basics and Extensions," From Fields to Strings: Circumnavigating 
Theoretical Physics, edited by M. Shifman, A. Vainshtein, and J. Wheater, (World Scientific, Singapore, 2005), Vol. I, pp. 
445-522, hep-th/0406216 

[4] G. V. Dunne and T. Hall, Phys. Rev. D 58, 105022 (1998). 



6 



[5] S. P. Kim, H. K. Loc, and Y. Yoon, Phys. Rev. D 78, 105013 (2008). 

[6] S. P. Kim and D. N. Page, Phys. Rev. D 65, 105002 (2002); Phys. Rev. D 73, 065020 (2006); Phys. Rev. D 75, 045013 
(2007). 

[7] G. V. Dunne and C. Schubert, Phys. Rev. D 72, 105004 (2005); G. V. Dunne, Q.-H. Wang, H. Gies, and C. Schubert, 

Phys. Rev. D 73, 065028 (2006). 

[8] H. Kleinert, R. Ruffini, and S.-S. Xuo, Phys. Rev. D 78, 025011 (2008). 

[9] W. Dittrich, Phys. Rev. D 19, 2385 (1979). 

[10] P. Elmfors, D. Persson, and B.-S. Skagerstam, Phys. Rev. Lett. 71, 480 (1993); P. Elmfors, P. Liljenberg, D. Persson, and 

B.-S. Skagerstam, Phys. Rev. D 51, 5885 (1995). 

[11] M. Loewe and J. C. Rojas, Phys. Rev. D 46, 2689 (1992). 

[12] A. K. Ganguly, P. K. Kaw, and J. C. Parikh, Phys. Rev. C 51, 2091 (1995). 

[13] J. Hallin and P. Liljenberg, Phys. Rev. D 52, 1150 (1995). 

[14] P. Elmfors and B.-S. Skagerstam, Phys. Lett. B 348, 141 (1995); Erratum, Phys. Lett. B 376, 330 (1996). 

[15] H. Gies, Phys. Rev. D 60, 105002 (1999); Phys. Rev. D 61, 085021 (2000). 

[16] S. P. Kim and H. K. Lee, Phys. Rev D 76, 125002 (2007). 

[17] C. M. Caves and B. L. Schurnakor, Phys. Rev. A 31, 3068 (1985). 

[18] B. L. Schumaker and C. M. Caves, Phys. Rev. A 31, 3093 (1985). 

[19] H.-Y. Fan and C. Li, J. Opt. B: Quantum Semiclass. Opt 6, S502 (2004). 

[20] N. B. Narozhnyi and A. L Nikishov, Yad. Fiz. 11, 1072 (1970) [Sov. J. Nucl. Phys. 11, 596 (1970)]. 

[21] S. P. Gavrilov and D. M. Gitman, Phys. Rev. D 53, 7162 (1996). 

[22] S. P. Gavrilov, D. M. Gitman, and J. L. Tomazelh, Nucl. Phys. B 795, 645 (2008). 



